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ABSTRACT 

Metallicity is expected to influence not only the lives of massive stars but also the outcome of 
their deaths as supernovae (SNe) and gamma-ray bursts (GRBs). However, there are surprisingly 
few direct measurements of the local metallicities of different flavors of core-collapse SNe. Here we 
present the largest existing set of host-galaxy spectra with H II region emission lines at the sites of 35 
stripped-envelope core-collapse SNe. We derive local oxygen abundances in a robust manner in order 
to constrain the SN Ib/c progenitor population. We obtain spectra at the SN sites, include SNe from 
targeted and untargeted surveys, and perform the abundance determinations using three different 
oxygen-abundance calibrations. The sites of SNe Ic (the demise of the most heavily stripped stars, 
having lost both H and He layers) are systematically more metal rich than those of SNe lb (arising 
from stars that retained their He layer) in all calibrations. A Kolmogorov-Smirnov test yields the very 
low probability of 1% that SN lb and SN Ic environment abundances, which are different on average 
by ^0.2 dex (in the Pettini & Pagel scale), are drawn from the same parent population. Broad- lined 
SNe Ic (without GRBs) occur at metallicities between those of SNe lb and SNe Ic. Lastly, we find 
that the host-galaxy central oxygen abundance is not a good indicator of the local SN metallicity; 
hence, large-scale SN surveys need to obtain local abundance measurements in order to quantify the 
impact of metallicity on stellar death. 
Subject headings: galaxies: abundances — supernovae: general 



1. INTRODUCTION 

Understanding the progenitors of the most ener- 
getic cosmic explosions, particularly gamma-ray bursts 
(GRBs) and supernovae (SNe), is a grand pursuit. 
Known already from their spectra is that stripped- 
envelope core-collapse SNe (CCSNe; "stripped SNe" 
hereafter) have progressively (from Types lib to lb 
to Ic) larger amounts of their outer hydrogen and 
helium envelopes removed prior to explosion (e.g., 
iClocchiatti et al.lll996t lFirippenkolll997l ). However, the 
dominant mechanism of stripping is not well known, 
nor are basic quantities such as the mass and metal- 
licity of their stellar progenitors. The exciting connec- 
tion between long-duration GRBs and broad-lined SNe Ic 
(SNe Ic-bl) and t he existence of SNe Ic-b l without ob- 
served GRBs (see iWooslev fc Bloom! 120061 for a review) 
raises the question of what distinguishes a GRB progen- 
itor from that of an ordinary SN Ic-bl without a GRB. 
Clear knowledge of the stellar progenitors of various ex- 
plosions is essential for understanding the endpoints of 
stars over a broad mass range and for map ping the chem- 
ical e nrichment history of the universe ([Nomoto et al.l 
[20061 . 

Two progenitor channels have been proposed for 
stripped SNe: either single massive Wolf-Rayet (WR) 
stars with main-sequence (MS) masses of > 30 M Q that 
have experienced mass loss during the MS and WR 
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stage s (e.g., iFilippcnko fc Sargentl [1985: Wo oslev et al.l 
1993), or binaries from lower- mass He stars that have 
been stripped of their outer en velopes through interac- 
tion ([Podsiadl owski et al.ll2004L and references therein), 
or a combination of both. Attempts to directly identify 
SN Ib/c progenitors i n pre-explosion images have not yet 
been successful (e.g ., iGal-Yam et al.l [20051 : iMaund et al.l 
l200ltlSina^l2009l) . 

A more indirect but very powerful approach is to study 
the environments of a large sample of CCSNe in or- 
der to discern systematic trends that characterize their 
stellar populations. SNe Ic tend to be found in the 
brightest regions of their host galaxies ([Kelly et al.ll2008|) 
and ar e more closely associated with H II regions than 
SNe II (|Anderson fc Jam es 2008, and references therein). 
SNe lb also tend to be found in bright regions of their 
respective hosts, less closely coupled than SNe Ic but 
more so than SNe II. This evidence suggests the pro- 
genitors of SNe Ib/c may thus be more massive than 
those of SNe II, which are ~ 8-16 M Q (see lSmarttll2009l 
for a review). Other studies attempt to measure the 
metallicit y by using the SN host-galaxy lumino sity a s 
a proxy ([Prantzos fc Boissier 2003; Arc avi et al.l 12010). 
or by using the metallicity of the galaxy center mea- 
sured from Sloan Digital Sky Survey (SDSS) spectra 
(iPrieto et al.ll2008f) to extrapolat e to that at the SN po- 
sition ([Boissier fc Prantzosl[2009[ ). 

Those prior metallicity studies do not directly probe 
the local environment of each SN (which is different 
from the galaxy center due to metallicity gradients), 
nor do they differentiate between the different SN sub- 
types. Here we present a statistically significant sample 
of stripped SNe (SNe lib and lb, SNe Ic, and SNe Ic-bl) 
with robust, uniform, and direct determinations of their 



local metallicity in order to quantify the impact of metal- 
licity on massive stellar de aths, building on our previous 
work (jModiaz et all 120081). We note that in the final 
stages of this research. lAnderson et al.1 (|2010( ) reported 
on a similar topic; we briefly note the differences in the 
sample and results below. 

2. THE SUPERNOVA SAMPLE AND ASSOCIATED 
HOST GALAXIES 

In Table [TJ we present the SN sample for which we 
measured local metallicities. It consists of 35 low-redshift 
(z < 0.18) stripped SNe, selected from the International 
Astronomical Union Circulars (IAUCs j3 according to the 
following criteria: (1) well-determined SN subtype, (2) 
discovered in targeted and untargeted surveys, and (3) 
H II region emission at the position of the SN (within 
our slit width of 1", see §131, as seen in our spectra, with 
which to directly determine the metallicity at the SN 
site. We discuss the potential impact of our selection ef- 
fects in § 15.11 While our sample is heterogeneous and not 
complete, we have good reasons to believe that it gives 
a fair representation of the kinds of environments that 
give rise to observed stripped CCSNe. Besides having 
SNe from tradi tional searches (e.g., t he Lick Observa- 
tory SN Search; iFilippenko et al.ll2001l) that target lumi- 
nous galaxies, we include SNe from untargeted surveys 
(e.g., SDSS, Nearby SN Factory; see IModiaz et all 120081 
for more detailed discussion) in order to mitigate any 
potential metallicity bias: since SN host galaxies in tar- 
geted searches are preferentiall y more luminous, the y are 
usually also more metal rich (jTremonti et al.ir2004f ). 

Furthermore, we include host-galaxy spectra of the 
broad-lined SNe Ic w ithout observed GRBs presented by 
iModiaz et al.l (|2008D . which were reduced and analyzed 
in the same fashion as the data presented here. For com- 
pleteness, we also include metallicity me asurements of 
all sp ectroscopically confirmed GRB-SNc (Modia z et all 
2008, and references therein; Christcnsen et all 2008: 



Chornock et "all l2010t [Starling et alJl2O10F 



The total sample of SNe lb, Ic, and Ic-bl whose local 
metallicities we are analyzing here amounts to 47 SNe 
(without observed GRBs). 

3. OPTICAL SPECTROSCOPIC OBSERVATIONS 

Optical long-slit (1" wide) spectra of the locations of 
faded SNe were obtained with the 10-m Keck I telescope 
using the Low R esolution Imaging Spectrometer (LRIS; 
lOke et all 119951 ) plus atmospheric dispersion corrector 
(ADC) on a number of nights 2007-2010. We generally 
employed a combination of the 300/5000 grism on the 
blue-side CCD and the 400/8500 grating on the red-side 
CCD Here we also include our Keck LRIS/ADC obser- 
vations of stripped CCSNe in those cases where super- 
imposed H II region emission lines were visible in the SN 
spectra (e.g.. IModiaz et al.ll2009l: iSilverman et al.H2009D . 

All optical spectra were reduced and calibrated with 
standard techniques in I RAF and our own IDL rou- 
tines for flux calibration (|Matheson et a l. 2008). For 
cases where the SN was still present, we eliminated 
the SN contribution following the successful method of 
IModiaz et al.l (|2008l ). After correcting all spectra for 

6 http://cfa-www.harvard.edu/iau/cbat.html . 



their recession velocities we measured optical emission- 
line fluxes by fitting Gaussians to the individual lines 
via the splot routine in IRAF. For the derivation of the 
statistical errors, which typicall y amount to 5-10% of 
the em ission-line flu xes, we foilowrPerez-Montero fc Dfazl 
(|2003! ) and, in part. lRupke et al.l (|2010t ). 

4. METALLICITY MEASUREMENTS 

The nebular oxygen abundance is the canonical choice 
of metallicity indicator for studies of the interstellar 
medium (ISM) , since oxygen is the most abundant metal, 
only weakly depleted, and exhibits very strong nebu- 
lar emission lines in the optical wavelength range (e.g., 
iTremonti et al.ll2004D . Using our measured line fluxes of 
[O II], [O III], [N II], Ha, and H/3, we correct for redden- 
ing via the Balmer decrement and the standard Ga lactic 
reddening law with Ry = 3.1 (jCardelli et al.l[l989T ). and 
compute the gas-phase oxygen abundance via strong-line 
diagnostics. 

We employ three independent and well-known cali- 
bratio ns: (1) the itera tive method of iKewlev fc Dopital 
(|200l ). as updated bv IKewlev fc Ellison! (l2008fl (hence- 
forth KD02-comb); (2) the calibration bv iMcGaughl 
(119911) (hencefo rt h M9 1); and (3) the diagnostic of 



iPettini fc Pagel (|200l (both PP04-O3N2 and N2) 
which is close to the direct electron temperature (T e ) 
scale. Moreover, we compute the uncertainties in the 
measured metallicities by explicitly including the sta- 
tistical uncertainties of the line-flux measurements and 
those in the derived SN host-galaxy reddening, and prop- 
agate them into the metallicity determination. Since 
the PP04-03N2 scale utilizes ratios of lines that are 
very close in wavelength, the effects of uncertain red- 
dening and scaling between the blue and red LRIS 
CCDs have negligible impact on the abundance measure- 
ments (something we tested). T he independently pub - 
lished metallicities of SNe 20 06jc (Pa storello et all2007|) , 
2007uy, and 2008D (jThone et al.ll2009D agree with our 
values within the uncertainties. 

5. RESULTS 

Figure Q] shows the cumulative distributions of local 
metallicities in each of the three scales for different types 
of stripped CCSNe: the ordinate indicates the fraction 
of the SN population with metallicities less than the ab- 
scissa value. The SN subtypes shown are SNe lb (includ- 
ing SNe lib) , Ic, and Ic-bl (without observed GRBs) , and 
their respective numbers for which the requisite emission 
lines for that diagnostic were available are given in the 
legend. Furthermore, we show confidence bands around 
each cumulative trend, which we computed via bootstrap 
with 10,000 realizations based on our metallicity mea- 
surements and their associated uncertainties. The metal- 
licity measurement s show well-known offset s between dif- 
ferent scales (e.g.. IKewlev fc Ellison! [2008h ; however, in 
each scale, SNe Ic are more likely than SNe lb to be 
found in metal-rich environments, while SNe Ic-bl (with- 
out observed GRBs) have environments that are similar 
in metallicity to those of both SNe lb and SNe Ic. A 
Kolmogorov-Smirnov (K-S) test reveals that the proba- 
bility that both the SN lb and SN Ic local host-galaxy 
metallicities have been drawn from the same parent pop- 
ulation is low: 1% (in PP04-O3N2), 7% (in KD02-comb), 
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Fig. 1. — Cumulative fraction (solid lines) of measured oxygen 
abundances at the SN position of different types of CCSNe with 
three different metallicity diagnostics and their confidence bands 
(dotted lines, see text). SNe Ic (the demise of the most heavily 
stripped stars that lost much, if not all, of both their H and He 
layers) are systematically in more metal-rich environments than 
SNe lb (SNe arising from less stripped stars that retained their He 
layer) . This is a robust conclusion since the trend is independent of 
the adopted metallicity diagnostic. The PP04-O3N2 scale is least 
impacted by reddening and flux-calibration uncertainties. Note 
that in this study, the SN lb subclass includes SNe lib as well. 



and 15% (in M91). The mean local metallicities in the 
PP04-O3N2 scale are 12 + log(0/H) = 8.49 ± 0.19 (for 
SNe lb), with a standard deviation of the mean (SDOM) 
of 0.012; 8.66 ± 0.12 with SDOM of 0.010 (for SNe Ic); 
and 8.47 ± 0.21 with SDOM of 0.016 (for SNe Ic-bl). In 
all scales, SN Ic sites have oxygen abundances that are 
0.14-0.20 dex higher than those of SN lb. 

Our results a re so mewhat at odds with those of 
lAnderson et al.l (|2010f) . who find a slight and insignif- 
icant difference between environmental metallicities of 
SNe lb and SNe Ic in their sample from only targeted 
SN searches. We have no SNe in common with their 
sample, and a detailed comparison of the two different 
samples is required to resolve the discrepancy, which is 
beyond the scope of this Letter. 



5.1. Tests for Possible Systematic Effects 

While there are selection effects that went into our 
heterogeneous sample, none of them is expected to af- 
fect SNe lb systematically more than SNe Ic, and hence 
should not be able to cause our observed trend. We 
checked that the SN survey mode does not explain the 
observed trend; indeed, both SNe lb and SNe Ic were 
drawn in almost equal proportions from both targeted 
and untargctcd surveys (Table 1), and the relative dif- 
ference in the metallicity of SN lb and SN Ic environ- 
ments is visible even when only comparing SNe from the 
same survey mode, albeit with more noise because of the 
smaller numbers of objects. Furthermore, we checked 
that both SNe lb and SNe Ic span comparable redshift 
ranges, with median redshifts of 0.015 (for SNe lb) and 
0.017 (for SNe Ic). However, the broad-lined SNe Ic in 
our sample extend to larger redshifts, with a median of 
0.043. 

While some surveys may have difficulty discovering 
SNe in the central cores of bright galaxies, this detection 
difficulty does not appear to affect one SN type more 
than another in our SN sample: from Table 1, most of 
the SNe included here were found far from the central 1", 
and those SNe that have offsets less than 1" comprise all 
types (SN lb, Ic, Ic-bl). The only strong selection effect 
in our sample is that we require H II region emission lines 
to be present at the SN position with which we can deter- 
mine the ISM oxygen abundance, meaning that the SN 
location has to have had a large amount of recent (i.e., a 
few million years) star-formation activity. However, this 
requirement affects SNe lb and SNe Ic equally since our 
objects sample the same redshift range. 

5.2. Supernova Progenitors with Metallicity- Driven 

Winds? 

A reasonable suggestion for why the environments 
of SNe Ic are more metal rich than those of SNe lb 
is that metallicity-driven w inds ()Vink fc de Koter] 120051 : 
ICrowther fc Hadficld 2006) in the progenitor stars prior 
to explosion are responsible for removing most, if not 
all, of the He layer whose spectroscopic nondetection 
distinguishes SNe Ic from SNe lb. This explanation 
may favor the single massive WR progenitor scenario 
as the dominant mechanism for producing SNe Ib/c 
(|Filippenko fc Sargent] [19851: IWooslev et all I1993D . at 
least for those in large star-forming regions (§ 15. ip . While 
the binary scenario has been sugg ested as the d ominant 
chann el for numerous re asons (see ISmartt 2009] for a re- 
view; iSmith et al.ll201ll ). we cannot ass ess it in detail, 
since none of the theoretical studies (e.g., Eldridgc et al. 
120081 and references therein) predict the metallicity de- 
pendence of the subtype of stripped SN. However, our 
result s are consistent with the suggestion of ISmith et all 
(|2011|) that SNe Ic may come from stars with higher 
metallicities (and masses) than SNe lb, even if they are 
in binaries. 

The distribution of the SNe Ic-bl is puzzling: SNe Ic- 
bl seem to occur at metallicities between those of 
SNe lb and SNe Ic, except GRB -SNe which are found 
at very low metallicit ies (but see iSoderberg et ail 120101 : 
iLevesaue et al.ll2010at also see lLevesaue et al.ll2010bl for 
other GRB host metallicities). This may indicate an- 
other key ingredient beyond metallicity for producing or- 



dinary SNe Ic-bl, perhaps magnetic fields. 

We find only a factor of 5 difference between the low- 
est metallicity Z for SNe lb and the highest metallic- 
ity for SNe Ic. Since the mass-loss ra te M is propor- 
tional to Z 086 (jVinkfc deKoterl 120051 ). this difference 
in Z would imply a maximum factor of 4 difference in 
M between SN lb and SN Ic progenitors. The question 
remains whether this small difference in M is enough to 
be responsible for removing all of the He layer, or other 
factors are responsible that have a higher dependence 
on metallicity than line-driven winds, or the metallic- 
ity trend simply c orrelates with another property such 
as pr ogenitor mass (|Kellv et al.l l2008; An derson fc Jame"sl 
2008) that may determine the SN outcome. All obs erva- 
tions and theoretical work (see lBastian et aLll2010l for a 
review) indicate that the initial mass function is univer- 
sal at the metallicities found herein. 

Our result s valid ate the independent hypothesis of 
lArcavi et al.l ( |2010f ). which was based on indirect data: 
as an explanation for the fact that none of the 15 CCSNe 
found by the Palomar Transient Factory (P TF) in dwarf 
galaxi es (Mr < — 18 mag) was a SN Ic, lArcavi et al.l 
(|2010l ) suggest that SNe Ic do not occur at low metallic- 
ity since low-luminosity galaxies usually have low metal- 
licity, in contrast to SNe Ic-bl of which two were found 
by PTF in dwarf hosts. Here we have supporting direct 
evidence; we show that SN Ic host environments have 
systematically higher metallicities than those of SNe lb, 
while those of SNe Ic-bl encompass both low and high 
abundances. Nevertheless, it is important to measure 
metallicities directly and not rely on the host-galaxy lu- 
minosity (L) as a proxy, as we show next. 

5.3. The Need for Local Metallicity Measurements 

Since we possess direct local metallicity measurements, 
here we test whether SN host-galaxy luminosity and mea- 
sured nuclear metallicity is a good proxy for the local 
SN metallicity, as assumed in some studies. To that end, 
we drew the SN host-galaxy luminosities from the Lyon- 
Meudon Extragalactic Database (HyperLEDAJJ as their 
sample has been homogeneously compiled, and adopt 
their reported absolute B-band magnitudes, Mb, in Ta- 
ble 1. Note that the lowest-luminosity SN host galaxies 
are from untargeted surveys. Furthermore, 14 galaxies 
in our sample have nuclear SDSS spectra and their PP04 
metallicities were taken from SDSffl 

In Figure [5] we plot the measured oxygen abundance 
(on the PP04-O3N2 scale) at the position of stripped 
CCSN vs. the SN host-galaxy luminosity, and for com- 
parison the oxygen abundance as would be inf erred from 
the SDSS L-Z relationship (jTremonti et al.1 120041 ) in- 
cluding lcr uncertainties. For co nsistency, we have con - 
verted the L-Z relationship from iTremonti et al.l (2004) 
to the scale of PP04-O3N 2 using the empirical calibra- 
tions of |KewIey^^jnhion| ((2008|) . Figure [2] demonstrates 
that nuclear metallicity when derived from the SN host 
luminosity is not a good proxy for the local oxygen abun- 
dance of the environments of SNe: the local metallicities 
are often lower than t he inferred central on es (because of 
metallicity gradients; Ivan Zee et al.l [1998). but also oc- 
casionally larger (e.g.. lYoung et al.1l2010f) . and in most 

7 http://leda.univ-lyonl.fr/ . 

8 http://www.mpa-garching.mpg.de/SDSS/DR7 



9.5 



CP 

o 

+ 

CM 



9.0 



c 
o 

'm 
o 
a. 



8.5 



o 8.0 



<D 
O 

c 
D 
"D 

c 

.Q 
< 

c 
a; 
en 
>^ 



7.5 



7.0 



-i — | — i — i — i — |- 



- 1 — i — i — 



- 1 — i — i — 



SN Ib + llb 
□ SN Ic 

SN Icbl 
XSN Icbl-GRB 



PP04-03N2 




* * 



* 



Systematic PP04-O3N2 error 



5DSS M B -Z Relationship in PP04-O3N2 — 

-J I I I I I I I I I I I I l_ 



1.0 



o 



0.1 



-16 -18 -20 -22 

Host Galaxy Mg [mag] 

Fig. 2. — Measured oxygen abundance at the position of different 
SN types vs. the SN host-galaxy luminosity, and for comparison 
the oxygen abundance as inferred from the SDSS L-Z relationship 
UTrem'onti et al.l l2004), converted to PP04-O3N2 (solid line) in- 
cluding la uncertainties (dashed lines). The nuclear metallicity as 
derived from the SN host luminosity using the SDSS relationship 
is not a good proxy for the local oxygen abundance of SN envi- 
rons. H ere, we adopt a solar oxygen abundance of 12 + log(0/H) 
= 8.69 ( Asplund et al. 2009), but note that the exact value of the 
solar abundance has no impact on our conclusions. SNe found in 
targeted SN surveys are designated by extra circles, squares, and 
diamonds. The representative systematic error of 0.14 dex for the 
PP04-O3N2 scale is shown in the bottom-right corner. 



cases more deviant than the la metallicity uncertainties 
(i.e., 0.16 dex) would indicate from SDSS. The differ- 
ences between the predicted central and the measured 
local metallicity values range from —0.4 dex to +0.5 dex 
(up to 3cr away from the SDSS L-Z value), with a median 
of 0.2 dex. Moreover, the nuclear PP04-O3N2 metallic- 
ities are deviant from the local values by an average of 
0.13 dex, with extremes as large as 0.24 dex, and always 
more deviant than our measured uncertainties. Thus, we 
conclude that nuclear metallicities are not a good mea- 
sure of local metallicities. 

6. CONCLUSIONS 

We present a statistically significant sample of 35 new 
host-galaxy spectra of stripped CCSN (SNe lib and lb, 
SNe Ic, and SNe Ic-bl), and in combination with pub- 
lished host-galaxy spectra, we perform robust, uniform, 
and direct determination of their local metallicity. The 
aim is to search for metallicity trends that may let us 
differentiate between various debated SN Ib/c progenitor 
scenarios. We find that the environments of SNe Ic are 
systematically more metal rich in three scales than those 
of SNe lb (on average by 0.14-0.20 dex in all scales), with 
a K-S test yielding very small probabilities that they are 
drawn from the same parent population. We also show 
that SNe Ic-bl (without GRBs) are intermediate to those 



other classes. 

For the future, we recommend this kind of detailed 
local metallicity study for all subtypes of CCSNe from 
galaxy-impartial, large, and deep photometric surveys 
(e.g., PTF, Pan-STARRS, Skymapper, LSST) in order 
to comprehensively understand the impact of metallicity 
on massive stellar death, as well as to compute cosmo- 
logically important parameters, such as SN rates, as a 
function of metallicity. 

Note added in proof. After our manuscript was ac- 
cepted, iLeloudas et al.l (|2011| ) appeared on the arXiv- 
server on a similar topic. While we corrected one SN 
identification based on their work, we are not able to 
make more comparisons at this time. 
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TABLE 1 
The Sample of Stripped-Envelope CCSNe 
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SN Type 
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Note. — SNc above the horizontal line: Local host-galaxy spectra arc presented and analyzed here for the first time. SNe below the horizontal line: we measi 
line fluxes and computed oxygen abundances from previously published spectra. 
a SN Discovery Type: T — SN host galaxy was targeted; Non-T — SN host galaxy was not targeted. 

Sample and data from Modjaz ct al. (2008), which uses the same technique as this work. 
c Remeasurcd from spectrum in lSahu et a.1.1 (2009). 

Measured from spectrum in Stritzingc r et al.| | |2009D . 



